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ABSTRACT 
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The effects of ocean surface waves on optical signal fluctuations and optical refraction 
are evaluated. The effects on signal fluctuations were inferred from the ocean surface wave 
spectral density and temperature fluctuation spectral density. The effects of ocean waves on 
optical refraction are manifested in the near-surface gradient of the temperature profile. A 
modified temperature profile incorporating a statistical mixing length scale is proposed. The 
research points to the importance of certain parameters that should be considered when 
optical signals skim the ocean surface. 
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L INTRODUCTION 



A. BACKGROUND 

The importance of environmental effects on communication links and radar systems 
has been recognized for the last few decades. Low altitude electromagnetic wave propagation 
over the ocean is strongly dependent on atmospheric conditions and sea state. It is essential 
to accurately characterize atmospheric turbulence in order to predict the propagation of 
electromagnetic signals. For a signal skimming the ocean surface, an understanding of ocean 
surface wave behavior should be given equal emphasis in electromagnetic wave propagation 
modeling. Specifically, it is necessary to relate the variation of the index of refraction to 
atmospheric turbulence and sea state. 

The refractive index at optical frequencies can be approximated by 



n-\ =77.6(1 +7.52x10* 3 A' 2 )(P/7)10' 6 , 



( 1 ) 



where P is the atmospheric pressure in millibars, T is the temperature in Kelvin and X is the 
wavelength of light in pm. To determine the variation in the index of refraction, it is sufficient 
to know the variation of temperature in the atmospheric surface layer. In the surface layer, 
the vertical profiles of wind speed, u, temperature, T and humidity, q, all have logarithmic 
forms £ =(u, T, q ), given by [Ref. 1] 



1 






( 2 ) 



where E,. is the flux scaling parameter, E, s is the value at the surface, z 0 ^ is the roughness 
length, k is the von Karman constant, ( = z/L, L the Monin Obukov length and rj/g the 
empirically derived stability correction functions. 

Equation (2) shows that the model results in large refractive index gradients in regions 
close to the sea surface. An experiment conducted by the Naval Research Laboratory in June 
1994 [Ref. 2] showed that the result of the large index gradient caused rays from the 
transmitter to be bent more than that required to establish optical continuity between the 
transmitter and the receiver. 

Unlike the atmospheric boundary layer above land, the boundary layer above the sea 
is continuously interacting with the ocean surface waves. The temperature profile of Equation 
(2) derived for an immovable surface may not be valid in the region of a few meters above the 
sea surface. Experimental data on signal heights that commensurate with ocean wave heights 
demonstrate deviations from the logarithmic distribution given by Equation (2) [Ref. 3], 

B. OBJECTIVES OF THE RESEARCH 

This research evaluates the effects of ocean surface waves on optical fluctuations and 
relates the optical refraction to the near-surface gradient of temperature and ocean surface 
waves. Measurements from an optical link are used to infer the temperature profile and hence 
the refractive index profile. These are used to assess the effect of the ocean waves. Based on 
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the experimental data, a proposed temperature profile is used in conjunction with a ray tracing 
program to validate optical continuity. 

C. ANALYSIS APPROACH 

The research will first examine the small-scale atmospheric fluctuations due to the 
presence of the ocean surface waves on the optical signal fluctuations. Computed ocean 
surface spectral density and air temperature fluctuation spectral density are used to establish 
the influence of the ocean surface waves. Having established the influence of the ocean waves, 
modification to the near-surface mean temperature gradients is proposed. This is done by 
incorporating a random scaling length to account for the random motion of the ocean waves. 
With the proposed temperature profile, a ray tracing code was used to establish optical 
continuity between the transmitter and the receiver. 

This report comprises two volumes. Volume I contains seven chapters covering the 
following details: The experimental set-up and instrumented flux buoy are described in 
Chapter H Chapter III presents the ray tracing code used in conjunction with the proposed 
temperature profile to validate optical continuity between the transmitter and the receiver. 
The analysis of ocean surface wave spectral density is presented in Chapter IV and the bulk 
aerodynamic method used to analyze the meteorological data is covered in Chapter V. The 
results of data analysis are presented in Chapter VI. Chapter VII summarizes the findings and 
recommendations of this research. The complete set of data plots is contained in Volume II. 
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H. EXPERIMENTAL SET-UP 



A. OPTICAL LINK 

In the experiment, an optical link at 0.83 pm was set up from Point Pinos at 
Lat./Long. (36° 38' 10.93 l’N, 121° 55' 45.444" W) to Marina State Beach at Lat./Long. (36° 
42 01.365" N, 121° 48’ 33.651" W) . The propagation path was approximately 12.8 km and 
the bearing of the transmitter from the receiver was 237°. A sketch of the optical path is 
shown in Figure 1 . 

The transmitter consists of a 0.83 pm GaAs laser source with an average power of 
35mW and intensity modulated at 4.2 kHz. It has a beam aperture of 50 mm and a beam 
divergence half angle of 1.3 mrad. The receiver has a 50 mm aperture, a laser line filter with 
a bandwidth of 40 nm centered at 0.83 pm followed by an Oriel Instrument Inc. battery 
operated silicon photodetector and pre-amplifier assembly. Block diagrams of the transmitter 
and receiver are shown in Figures 2 and 3 respectively. During the experiment, the heights 
of the transmitter and receiver were varied and the received signal recorded on tape. 

The recorded signal was demodulated with a Stanford Research Systems lock-in 
amplifier with sensitivity set at 50 mV, time constant of 1 ms and internal oscillator set at 4.2 
kHz. The output of the lock-in amplifier was connected to a Labview National Instruments 
12 bit analog to digital card plugged into a Macintosh computer. The analog signal was 
sampled at 400 Hz. A total of 41 signal samples covering the period from 30 May 96 to 5 Jun 
96 were taken. 
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Figure 2. Transmitter Block Diagram 



D. 
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Figure 3. Receiver Block Diagram 



B. FLUX BUOY 



The flux buoy was deployed at approximately the mid point of the optical link, at 
Lat./Long. (36° 40 N, 121° 52’ W) where the ocean depth is about 69.5 m. The flux buoy was 
deployed to record meteorological and wave data. The flux buoy system contains two 
separate sensor suites, viz., a mean meteorology package and a high frequency sampling ‘flux’ 
package. Two sets of data were collected; the high frequency data (flux) of air temperature 
and wind velocity fluctuations sampled at 21 Hz, wave motion data sampled at 10.5 Hz and 
low frequency data (mean meteorology) of wind speed, wind directions, humidity, pressure, 
air and sea temperature sampled at one minute interval. A picture of the buoy is shown in 



The high frequency air temperature and velocity were measured by the Solent 
ultrasonic anemometer at 5.09m above mean tide level. The measurement technique is based 
on sound pulses transmitted back and forth between transducers. The time of flight in each 
direction, say t , and t 2 are measured. If c is the speed of sound, x the distance between the 
transducers, and there is an air flow v along the line of the transducers, the following 
relationships are readily derived: 



Figure 4. 




( 3 ) 



x 

(c-v) ’ 



( 4 ) 
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v = 0.5 x 



(5) 



1 1 



and 



c = 0.5 x 




( 6 ) 



By arranging three pairs of transducers in different orientations, the direction and magnitude 
of the incident air flow may be unambiguously derived. The air temperature is calculated from 
the speed of sound as follows [Ref. 4]: 



c 2 = yRT(l +0.51?) , 



(7) 



where y R is 403 m 2 s' 2 K' 1 for air, T is the air temperature in degrees Kelvin and q is the 
specific humidity in kg of water vapor per kg of air. 

The ocean surface wave data was measured by the Systron Donner Motion Pak. The 
Motion Pak consists of three accelerometers measuring the accelerations of the x, y and z 
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directions (sway, surge and heave) and three rate gyros measuring the angular rates, 0, § and 
i|r (pitch, roll and yaw). The data was sampled at 10.5 Hz. The Motion Pak uses a right- 
handed coordinate system for the angular rates about the x, y, and z axis. It was mounted so 
that the x-axis points towards the port side of the buoy, the y-axis points towards the bow and 
the z-axis downward. The z-axis voltage was adjusted so that when the buoy is at rest on a 
flat surface, the z-axis linear acceleration voltage output read zero voltage. 

Three wave staffs were deployed to measure the high frequency ocean surface waves. 
The three wave staffs were arranged such that each staff was located at the vertex of an 
equilateral triangle. The arrangement allows for the computation of the wave directional 
spectra. Unfortunately, wave staff number 2 was not functioning properly during the 
experiment. 

The buoy’s heading relative to magnetic North can be determined from the three 
Syntron compass outputs. The magnetic deviation from the true North at the buoy was 15° 
15 East. The buoy’s heading relative to the true North can be determined by subtracting the 
deviation from the buoy magnetic heading. 

The low frequency wind speed was measured by the R.M. Young Wind Monitor 
located at 3.80 m above mean tide level. The wind speed sensor is a helicoid shape propeller 
molded of polypropylene plastic. Threshold sensitivity of the propeller is 1.0 m/s. Rotation 
of the propeller produces a linear DC voltage signal directly proportional to the wind speed. 

The Everest IR sensor was used to measure the sea surface temperature. The IR 
radiation from the sea surface is collected by an optical lens system and focussed on an IR 
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detector. The detector in turn converts the IR signal to an electrical signal proportional to the 
sea surface temperature. 

The low frequency temperature and humidity were measured by the Rhotronics 
relative humidity and temperature sensor. The humidity sensor is a small hygroscopic 
capacitor that changes its value as a function of both water vapor pressure and temperature. 
The temperature sensor uses a platinum 100 ohm resistance temperature device (RTD) to 
measure air temperature. 

The flux data were stored in an Onset Model 6 Toshiba computer with 500 MB of 
disk storage on board the buoy. Each word is a two-byte data stored as ‘big endian’ 
hexadecimals (most significant byte first). The mean meteorology data were stored in another 
Onset Model 6 computer with 200 MB of disk storage. The data were stored in ASCII 
format. The bulk aerodynamical method was used to process the mean meteorological data. 
C. ENVIRONMENTAL INFLUENCE AND LIMITATIONS 

The Monterey Bay has deep water with ocean-like swells. The optical link from Point 
Pinos to Marina State Beach crosses Monterey Bay where there is relatively very little 
shoreline influence. The buoy was located at the mid point of the link in the deepest water and 
farthest from land with open-ocean wind and waves. 

The limitation in the deployment of the buoy at a single location leads to the 
assumption of a spherically stratified atmosphere with a lateral homogeneous refractivity 
along the path of propagation. This assumption is reasonable only in a well-mixed 
environment over the ocean. 
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In the experiment, the spatial average optical measurement is compared with the 
temporal average meteorological data. In a situation where there is lateral homogeneity along 
the mean wind direction, the two averages could be considered to be equivalent [Ref. 1], 
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Figure 4. Flux Buoy 
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m. RAY TRACING 



The ray tracing code was used in conjunction with the standard logarithmic 
temperature profile and the proposed modified temperature profile to trace rays emanating 
from the transmitter to the receiver. Optical continuity is said to be achieved if there exists 
at least one ray above the ocean wave height connecting the transmitter to the receiver. 

The ray tracing technique used here is based on the methodology developed by Lehn 
[Ref. 5], The model assumes that the atmosphere consists of a number of concentric spherical 
layers where the temperature is assumed to vary linearly within each layer. In this case, the 
curvature of a nearly horizontal ray is approximately constant within each layer. Replacement 
of circular arc by parabolic arc leads to a simpler quadratic expression for the ray trace. 

For a ray traversing a stratified atmosphere above a flat earth, Snell’s Law dictates 

that 

nsinO = A, (8) 



where n is the refractive index of the air at any point on the ray, 0 is the angle between the 
ray and the vertical, and A is a constant for the ray. 

In rectangular coordinates, the curvature, C, is given by 
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C = ± 



( 9 ) 



d 2 z 



dx 7 



dx 



3/2 



where x is the horizontal distance and z is the elevation above the earth’s surface. The slope 
of the ray is given by 



dz _ cos0 
dx sin0 

( W 2 — ^ 2 ) 1/2 ( 10 ) 

A 



and the second derivative is given by 



dh 

dx 2 



d_ 

dz 



dz 

dx 



dz_ 

dx 



n dn 
A 2 dz 



( 11 ) 



With Equations (10) and (11), the ray curvature becomes 



C = - 



sin0 dn 
n dz 



( 12 ) 
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The sign in Equation (12) has been chosen so that a positive curvature corresponds to a ray 
that is concave towards the earth. At optical frequencies, the refractive index of air can be 
expressed as 



n = 1 + 



e^£ 

T ’ 



(13) 



where p is the pressure, T is the absolute temperature in degrees Kelvin, e = 226x1 O' 6 m 3 kg' 1 
and p = 3.48xl0' 3 J' 1 kg K. Note that Equation (13) is equivalent to Equation (1) when the 
dependence on wavelength is set at X = 0.83 pm and the pressure expressed in N/m 2 . The 
derivative of n is given by 



£ £ = _ €gp Q dT 
dz T 2 dz ’ 



(14) 



where the relation p = pp/Thas been used, p is the density of the air and g is the acceleration 
due to gravity. Substituting Equation (14) into the expression for the ray curvature leads to 



c = sine €pp L 1 dT 
n T 2 [ dz 

= e P 

(1+ep) 




(15) 
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In ray tracing, the circular arcs that represent the earth’s surface, the ray and the layer 
boundaries are all approximated by parabolic arcs. The coordinate system is shown in Figure 
5. The surface of the earth is given by 



z = 




(16) 



where R is the radius of the earth. Similarly a layer located at a height h above the earth’s 
surface is given by 



z 



— +h 

2(R+h) 




(17) 



An optical ray with a radius curvature, r = 1/C, starting at x = 0, z = z 0 and making an angle 
of 4> 0 with the horizontal is given by 



z 




+ xtanc{> 0 + z 0 . 



(18) 
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Within the f layer, the ray path is given by Equation (18) with a radius of curvature 
r, . Ray tracing is performed by determining the intersection of the ray with the layer 
boundaries. The i h layer lies between elevations h imI and A , . A ray can enter a layer either from 
above or below and could exit the layer either from top or bottom of the layer. Figure 6 
shows a ray entering the i lh layer from lower boundary. It will proceed to P if r i > R or r, < 
0 and to Q if r is positive and sufficiently small. The coordinates of P are found by solving 
Equations (17) and (18) with A = A,. in Equation (17) and z 0 = A,. ; in Equation (18). The x- 
coordinate of P is 



If the discriminant is negative, the intersection P does not exist. Then Q has to be determined. 
This is done by letting A = A,., in Equation (17). The x-coordinate of Q is 



Similarly, for a downward-heading ray, the x-coordinate of P at layer boundary A,. ; is given 




(19) 



27?/' J tan<t> 0 



R-r : 



( 20 ) 
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by 




/ 

-tan(j> 0 - tan 2 <j> 0 - 



2(h-h i _ l )(r-R)y / > 



( 21 ) 






Again, a negative discriminant implies that P does not exist and the ray bends upwards and 
returns to layer boundary h { at Q. The x-coordinate of Q is determined from Equation (20). 

New coordinates are used each time the ray crosses a layer boundary. Figure 7 shows 
the case of an upward-heading ray. As the curvature is small, the x-coordinate of the new 
origin is taken to be u P . The layer boundaries can still be described by Equation (17). The 
elevation angle of the ray is given by 



where m r and m [ are the ray and layer slopes respectively at P. The slopes are found by 
differentiating Equations (17) and (18) giving 



4> 0 = tan“‘w r - tan' 1 /^ , 



( 22 ) 





( 23 ) 
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Figure 5. Ray Tracing Coordinate System [From Ref. 5] 
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Figure 6. Upward-heading Ray in Layer / [From Ref. 5] 
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Figure 7. Change of Coordinates when the Ray enters a New Layer [From Ref. 5] 
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IV. OCEAN WAVE SPECTRAL ANALYSIS 



A. MEASUREMENT APPROACH 

In order to establish the possible effects of ocean surface waves on the optical link, 
it is necessary to examine the ocean surface wave spectra and wave heights at the time when 
each sample of optical signal was measured. The wave information was collected by the 
motion pack and wave staff. The motion pack consists of three accelerometers measuring the 
accelerations of the x, y and z directions (sway, surge and heave) and three rate gyros 
measuring the angular rates, 0, <f> and i|/ (pitch, roll and yaw). The data was sampled at 10.5 
Hz. 

The true surface displacement is obtained by summing the following three 
components, viz., the wave height as measured by the wave staff, z m , the displacement of the 
wave staff due to the translational motion of the buoy, z , , and the displacement of the wave 
staff due to the rotational motion between the wave staff and the motion pack [Ref. 6], 

As the wave staff and the linear accelerometers were mounted on the moving 
reference of the buoy, it is necessary to compute the surface displacement in an earth- 
referenced system. The conversion of the buoy-referenced system to the earth-referenced 
system can be achieved via a transformation matrix, T be [Ref. 7] given by 
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BE 



cos0 cosi|j sin<(> sin0 cosxjr - cos<() simjx 
cos0 simjf sincj) sin0 sini}/ + cos<J> cost}; 
-sin0 sin())cos0 



cos(j) sin0 cosi}/ + sine}) simjr 
cose}) sin0 sim}r - sine}) cosi|r 
cosc{)cos0 



• (24) 



If a B represents the translational accelerations measured by the strap-down 
accelerometers on the buoy, the earth-referenced accelerations, a E are given by 



a E T be a B + & ’ 



(25) 



where g = ( 0 , 0, -g) is the acceleration due to gravity. The earth-referenced vertical 
displacement due to the translational motion of the buoy is given by 



= JJl a B ’ ( - sin0,cos0sin(}),cos0cos(}))-g]<#dif 



(26) 



Similarly, the earth-referenced vertical displacement measured by the wave staff is given by 



Z mE = ^COS0COS<}) . 



(27) 
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Finally the displacement due to the rotational motion is determined by integrating the 
rotational velocity vector, Q x Tbe Lb where Q is the angular velocity of the buoy measured 
by the rate gyro and Lg is the distance between the accelerometer and the water surface at the 
wave staff. The angular velocity, Q is given by 



Q = 



-0 sim|r + <j>cos0 costjf 
0 cost); + <j> cos0 sintjt 
ijr-<f)sin0 



(28) 



The vertical displacement due to the rotational motion is given by 



z rE = f [L E2 ( -0sim|r + 4> cos0 cosij; ) - ( 0 cosijr +4>cos0sinil;)]^/ , (29) 



where 



^be L b 



(30) 



and the subscripts El and E2 represent the component of the vector. 
The true ocean surface displacement is given by 
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( 31 ) 




Once the ocean surface displacement is determined, the rms ocean surface displacement, o 
can be computed and the significant wave height, H s is given by 



B. ANALYSIS OF OCEAN SURFACE DISPLACEMENT 

The analysis of ocean surface displacement or elevation was carried in the frequency 
domain using Fourier Transform. The rate gyros do not have good low frequency response 
and therefore it is necessary to use complementary filtering to patch the low frequency tilt 
reference of the rate gyros with the linear accelerometers as follows 



H = 4o . 

S 



(32) 



0 = Q Hp -sm , 



(33) 



4 > = %P + sin ' l ( a x / S) L p 



(34) 



and 
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(35) 



i|/ = y]f Hp + sinr\com) Lp , 



where a# and a y represent the observed linear acceleration in the x and y direction, the 
subscripts HP and LP represent the high pass and low pass filtered output and com is the 
compass output. The first term on the right hand side of Equations (33) to (35) represents the 
high pass angular motions from the rate gyros which provides the wave-induced angular 
motions and the second term represents the low frequency tilt reference from the 
accelerometers. The effect is an all-pass filter which removes the unwanted drift in the rate 
gyros while retaining the low frequency tilt reference. 

To compute z rE and z# , the integration and double integration are done in the 
frequency domain by dividing the Fourier Transform of z rE and z $ with (ico) and (ioo) 2 
respectively. The division by o> introduces a singularity at to = 0. However, there is no ocean 
wave with frequency lower than 0.04 Hz, i.e., 25 second waves, since there is not enough 
fetch in the ocean to sustain the waves without breaking. The waves in the Atlantic Ocean and 
Pacific Ocean rarely exceed 20 seconds and 25 seconds respectively [Ref. 8], The singularity 
problem is alleviated by performing a high-pass filtering of the wave spectrum at 0.04 Hz. The 
power spectral density of the ocean surface waves and significant wave height were computed 
for the duration of time when the optic signal samples were being measured. 
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V. BULK AERODYNAMICAL METHOD 



The bulk aerodynamical formulations for the fluxes of momentum, sensible heat and 
latent heat across the air-sea interface and vertical profile estimation are given by 

t = P C D (z)[u(z)-u$ , ( 36 ) 



Qh = -9 c p C^z)[u{z)-u^[T(z)-T^ (37) 



and 

Q e = ~? L e C E {z)[u(z)-u^[q{z)-q} , ( 38 ) 

where p is the air density, c p is the specific heat of air, L e is the latent heat of vaporization, C D , 
C H and C E are the height dependent bulk transfer coefficients for momentum, sensible heat 
and latent heat respectively, u, T, q are the wind speed, potential temperature and specific 
humidity respectively, z is the height above the sea surface and the subscript s denotes the 
respective values at the sea surface. 

In the atmospheric surface these fluxes can be expressed according to Monin-Obukov 
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Similarity (MOS) theory scaling parameters as 



T 




( 39 ) 



2 * 



-? c p u J \ 



(40) 



and 



Q e = ~ pL e u m q , , 



(41) 



where T, and q* are the MOS surface layer scaling parameters for momentum, temperature 
and humidity respectively. By combining Equations (36) to (41), these scaling parameters can 
be defined according to the bulk transfer coefficients as 



«. = C^ 2 [u(z)-u s ] 



(42) 



T. = C„Cv m {T(z)-T,} 



(43) 
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and 



?. = C D m [q(z)-q,] 



( 44 ) 



The bulk transfer coefficients can be defined as 



C D iz) = 



[ln(z/zJ-t|r a (z/L)] 2 



( 45 ) 



CtfOO = 



^ 

[In (z/zj -^ z/L)]Mz/z 0T ) -i^z/I)] 



( 46 ) 



and 



C e( z ) = 



K 2 

[ln(z/z 0K ) - iJ;(z/L)][ln(z/z 0g ) - y q (z/L)] ’ 



( 47 ) 



where k is the von Karman constant, and t|r a , i|r r and are the integrated forms of the 
dimensionless stability correction functions for wind speed, temperature and specific humidity 
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respectively, given by 



V U WL) = 



21n[(l +x)/2] +ln[(l +x 2 )/2] - 2tan _1 (x) + ti/2, z/L < 0 
-8 z/L, z!L > 0 ’ 



(48) 



t|r r (r/Z) = i|r 9 




21n[(l +x)/2], 
-lz!L, 



z/L<0 
z/L > 0 ’ 



(49) 



x = (1-20 z/L) VA 



(50) 



and z 0u , z 0T and z 0q are the roughness lengths for momentum, temperature and humidity 
respectively. The Obukov length scale, L is defined as 



L = 



[Kg(7>0.6ir<7,)J ’ 



(51) 



where g is the acceleration due to gravity and T v is the virtual temperature introduced to take 
into account the effects of humidity: 
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T v = 7X1+0.61?) 



(52) 



To determine the vertical profile of wind speed, temperature and specific humidity, 
iterative computation of u*, T., q,, z^, z 0T , z 0q and L are performed and the vertical profile 
calculated as 






u = u s + — [\n{zlzJ-1? u {zlL)} , 

1C 



(53) 



T=T+ — [ln(z/z or ) - tyrizIL)] 



K 



(54) 



and 



[In (zlz„ q )-^ q (zlD] . 

1C 



(55) 
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The temperature profile estimated in this manner was used as input to the ray tracing 
program to determine if optical continuity between the transmitter and the receiver can be 
established. 
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VL RESULTS OF DATA ANALYSIS 



A. OCEAN SURFACE WAVE SPECTRA AND WAVE HEIGHTS 

The power spectral densities of the ocean surface waves were computed for the 
duration when the optical signals were being measured. A total of 41 samples were computed. 
A typical sample of the wave spectrum is shown in Figure 8. It can be seen from the spectrum 
that the ocean surface waves are dominated by swells with a peak frequency of about 0. 1 Hz, 
i.e. 10 second waves. There are wind generated waves with a peak frequency at about 0.4 Hz, 
i.e., 2.5 second waves. The spectrum exhibits a high frequency roll-off as o' 4 . This is the 
established behavior of the equilibrium range of wind-sea spectra [Ref. 9 ]. 

The rms wave height, a, was computed from the power spectral density by the 
following relation 



N 2 /=i 



(56) 



where Z(j) is the Fast Fourier Transform (FFT) of the ocean surface displacement, z(t) and 
N is the sample size. For the 41 samples, the rms wave height ranges from 0. 18 m to 0.7 m. 
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Figure 8. Ocean Surface Elevation Spectrum. The straight line represents the -4 slope. 
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B. 



TEMPERATURE SPECTRAL DENSITY 



The aspect of atmospheric turbulence of relevance to optical propagation is the 
variation of the index of refraction or equivalently, the variation of the density of air. Air 
density depends on pressure and temperature. Pressure fluctuations are rapidly dissipated and 
their effects are negligible compared to temperature fluctuations. Thus, for optical 
propagation purposes, temperature fluctuations are used to characterize the optical effects 
of turbulence. 

In the experiment, temperature fluctuations sampled at 21 Hz were measured by the 
sonic anemometer at 5.09 meters above mean sea level. A typical plot of the power spectral 
density of temperature fluctuations is shown in Figure 9. The spectrum shows the classical 
-5/3 slope at high frequencies. This dynamic condition in the high frequency range agrees 
satisfactorily with theoretical predictions and measurements above land [Ref. 10], At lower 
frequencies, the spectrum shows evidence of ocean wave influence. Figure 10 shows the 
temperature spectral density and corresponding ocean surface spectral density of the wind 
waves. The spectral density of wind waves was obtained by performing high pass filtering on 
the swell-dominated wave spectrum, thus removing the swells from the spectrum. It can be 
seen that the temperature spectrum displays the characteristic spectral density peak at the 
peak of the corresponding wind wave spectrum indicating the influence of the ocean waves. 
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Figure 10. Temperature Spectrum and the Corresponding Wind Wave Spectrum. 
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c. 



OPTICAL SIGNAL FLUCTUATIONS 



In the analysis of the optical signals, the signal fluctuations are described by the log- 
normal variance although the variance of the signals can be used as a measure of fluctuations. 
The research does not attempt to relate the log-normal variance to the saturation 
phenomenon. The objective is to analyze how the effects of ocean surface waves affect the 
signal fluctuations per se. 

1. Log-normal Distributions 

To observe the signal fluctuations, a histogram plot of the normalized log intensity 
(normalized to the mean intensity) was performed for each signal sample. From the histogram 
plot, it was found that for signal samples where no apparent ocean wave occultation is 
observed, the intensity fluctuations exhibit a log-normal distribution. For signal samples where 
wave occultation is evident, deviation from the log-normal distribution results. Typical 
histogram plots of samples that exhibit log-normal distribution and ‘non log-normal 
distribution’ are shown in Figures 11 and 12 respectively. 

Log-normal distributed signals are known to exist in signals propagating through 
random inhomogeneous medium [Ref. 11], The distribution characteristics are associated with 
atmospheric radio noise or atmospherics as it is sometimes called. It has been found 
empirically [ Ref. 12] that the distribution P(E>R), where E is the instantaneous value of HF 
atmospheric noise and R is threshold voltage at the receiving point, approaches a Rayleigh 
distribution for small amplitudes. This is due to the interferences of many mutually 
overlapping small atmospherics with uniformly distributed phase. However for large 
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amplitudes, the distribution deviates considerably from the Rayleigh distribution and 
approaches a log-normal distribution. This is believed to be due to strong atmospherics that 
do not overlap in time. 

The theory of the distribution of peak value of atmospherics can be explained as 
follows. The peak value of atmospherics, E p can be represented by the following expression 





(57) 



where P 0 is the peak power of the atmospheric discharge, d is the propagation path, d j is a 
section of the propagation path, 3j is the attenuation constant of the j th section and AT is a 
constant. By writing Equation (57) as 



E p = e A = exp 



( n n \ 

-E V/ + ^ ln/> 0 - l»E<*y - In* 



/= 1 



/= 1 



(58) 



the sum of the attenuations of each section, A, is defined. It is assumed that the fluctuations 
of P 0 and d are negligible compared to the first term in the exponential of Equation (58). 
Thus, only the first term will thus be used to determine the distribution of A. Since most 
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atmospherics are propagated from long distances away, it is reasonable to assume that the 
attenuations 5^ are independent and none of the 8 y d? s make a predominant contribution 
to the sum of the attenuations. Therefore A is normally distributed with probability density 
function given by 



P( A) = 




exp 



(A-p) 2 

2o 2 



( 59 ) 



2 • 

where p is the mean, a A is the variance. From Equation (58) we have 



A = In E p 



( 60 ) 



and substituting Equation (59) into Equation (60), the density function of E p is 



KE p ) = 




0^-p) 2 



2o. 



( 61 ) 



Equation (61) is the log-normal distribution which has been found experimentally to apply to 
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the propagation of optical signals through atmospheric turbulence in the absence of 
occultation of the optical signals by ocean surface waves. 
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Figure 11. Optical Signal Intensity and Distribution Plots and the Corresponding Ocean Waves. 
‘...’ in ocean waves plot indicates height of optical ray path above the waves. 
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Figure 12. Optical Intensity and Distribution Plots and the Corresponding Ocean Waves Plot 
(with Occultation of Optical Signal by Ocean Waves). 
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2. Log-normal Variance and Mean 



From the above equations, the relationship between the log-normal mean intensity and 
log-normal variance can be derived. In general, an electromagnetic wave field E(r,t) can be 
represented by 

E{r,t) = I(j) m exp [ik-r - rot] , (62) 



where co is the temporal frequency of the electromagnetic field, k is the propagation vector 
and I(r) is the irradiance or intensity. In the presence of turbulence, Equation (62) becomes 

E(r,t) = <I(r)> V2 exp [ik-r - iut + x(fj) + <j>(r,0] > ( 63 ) 



where <I(r)> m is the average non fluctuating part of the irradiance, x(r,t) describes the 
fluctuating part and 4>(r, t) is the phase variation. The random variables, %(r, t) and 4 >(r, t) 
have Gaussian probability distributions. To examine the irradiance fluctuations, the 
instantaneous irradiance can be written as 



I(r,t) = <7(r)>exp[2xM] . 



(64) 
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It follows that I(r, t) has a log-normal distribution. Taking the average on both sides of 
Equation (64), and using the relation valid for any real-valued Gaussian random variable, 



( e az ) = exp 



a(z) + —a 2 o 2 

2 



(65) 



Equation (64) becomes 



</> = </>exp(2<x> + 2 a 2 ) . 



( 66 ) 



In the presence of turbulence, </> should be a constant independent of the turbulence. It 
follows that 



<X> = 




(67) 



A plot of log normal mean intensity against the log normal variance for those signal 
samples not occluded by the ocean surface waves in Figure 13 further verifies the log-normal 
distribution characteristics of the optical signals. 
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Figure 13. Log-normal Variance vs. Log-normal Mean. 
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D. EFFECTS OF OCEAN SURFACE WAVES ON OPTICAL SIGNAL 
FLUCTUATIONS 

Due to the random motion of the ocean waves, it is necessary to establish a way to 
measure the distance of the ray path above the ocean waves. The method adopted here is to 
calculate the normalized mean square distance, h 2 , (normalized to the variance of the ocean 
surface displacement) of the ray path above the ocean waves, given by 



h2 = -7 E (hr 2 ) 2 



( 68 ) 



where o 2 is the variance of the ocean surface displacement, z hl is the lowest height of the 
straight line ray path above mean sea level, z, is the instantaneous ocean surface displacement 
and the summation is performed over the duration of the optical signal transmission. A large 
mean square distance implies a lower probability of the ray path being blocked by the ocean 
surface waves. The normalized mean square distance is plotted against the log-normal 
variance of the signals. For sufficiently large mean square distance, i.e., for the mean square 
distance greater than 35, the variance of the optical signal increases as the mean square 
distance decreases. Those signal samples with lower mean square distance also display the 
same trend but with increased variability, i.e., the signal sample points are widely scattered. 
A plot of log-normal variance of the optical signal against hi 2 is shown in Figure 14. The signal 
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samples are numbered from 3 to 43 for easy identification of any signal sample. The list of 
optical signal samples is contained in Appendix A. 

From Figure 14, it can be seen that there is an apparent ‘peak’ in the signal 
fluctuations for signal samples 14 to 17. The higher signal fluctuations are due to the 
increased roughness of the ocean wind waves associated with the signal samples. The wave 
age, defined as the ratio of the ocean surface waves phase speed, C p to the wind speed, u is 
a good measure of roughness of the surface waves. Donelan, et. al. [Ref 13], demonstrated 
that younger waves traveling slower than the wind are rougher than mature waves. A plot of 
the optical signal log-normal variance against wave age for signal samples 3 to 5 and 14 to 
17 is shown in Figure 15. The plot indicates a decrease in signal fluctuations with increased 
wave age. However, signal samples 14 and 15 have higher signal fluctuations than samples 
3 and 4 although they are associated with wind generated waves of approximately the same 
wave age. The higher fluctuations of samples 14 and 15 are due to the influence of the swell. 
The influence of the swell on sample 14 and 1 5 is manifested in the temperature fluctuations 
spectral densities where the swell frequencies modulate the temperature spectra. This effect 
is shown in Figures 16 and 17. The corresponding temperature spectra for samples 3 and 4 
do not display any prominent peak in the swell frequencies. The temperature spectra and the 
ocean surface elevation spectra for samples 3 and 4 are shown in Figures 18 and 19. The 
presence of swell caused a change in the spectral distribution of the various components of 
the high frequency wind waves. There is experimental evidence that this leads to an increased 
in roughness of the ocean surface [Ref. 14]. 
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Two alternative analyses were performed to verify the trend of the signal fluctuations. 
In the first method, it is assumed that the ocean surface elevation follows a Gaussian 
distribution, T|(0,a). This is true for a sample with a large number of.points. Based on the 
height of the ray path above mean sea level, the probability of ocean surface waves protruding 
into the ray path can be computed from 



The percentage of the ocean waves protruding into the ray path is found by multiplying 
Equation (69) by 100. The percentage values are plotted against the signal log-normal 
variance. It was found that the signal fluctuations exhibit the same trend. A plot of the 
percentage of ocean waves protruding into the optical ray path is shown in Figure 20. 

In the second method, the log-normal variance of the signal samples are plotted 
against the difference between the ray path height and one significant wave height above the 
wave crests. The same trend in signal fluctuations was observed. For signal samples that are 
less than one significant wave height above the ocean surface waves, the sample points are 
widely scattered as in the above two analyses. The plot is shown in Figure 21. The analyses 
led to the conclusion that the random motion of the ocean waves is affecting air turbulence 




( 69 ) 
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immediately above the ocean surface and causing increased variability in signal fluctuations 



as the ray path gets closer to the ocean surface. 
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Figure 15. Optical Signal Log-normal Variance vs. Wind Wave Age. 
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Figure 16. Optical Signal Sample No. 14, Temperature Spectrum and the 
corresponding Ocean Surface Wave Displacement Spectrum 
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Ocean surface elevation spectrum, 31 May 96, 1115 hr 
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Figure 17. Optical Signal Sample No. 15, Temperature Spectrum and the 
Corresponding Ocean Surface Wave Displacement Spectrum. 
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corresponding Ocean Surface Wave Displacement Spectrum 
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Figure 19. Optical Signal Sample No. 4, Temperature Spectrum and the 
corresponding Ocean Surface Wave Displacement Spectrum 
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Log-normal variance 



Log-normal variance vs. percentage of ocean waves protruding into ray path 




Figure 20. Percentage of Ocean Waves Protruding into Optical Ray Path. 
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Log-normal variance 



Log-normal variance vs. ray path 1 sig wave height above ocean waves 




above Ocean Waves. 
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E. OPTICAL SIGNAL MODULATION BY OCEAN SURFACE WAVES 

The intensity spectrum for each optical signal sample was computed and compared 
with the corresponding ocean surface wave spectrum. The plot (Figure 22) shows that for 
optical signals that were partially occluded by the ocean surface waves, the optical signal 
spectrum displays peaks at the swell frequency of the ocean waves. The optical signal is said 
to be modulated by the ocean waves. 
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Figure 22. Optical Signal Spectrum and the Corresponding Ocean Wave Spectrum. 
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F. 



NEAR-SURFACE MEAN GRADIENTS 



In the above, only small-scale variations are considered whereby turbulence produces 
random fluctuations in the density of air and hence fluctuations in the refractive index. On 
larger scale variations, the main concern is on phenomena that affect the average position of 
the optical beam as it traverses the atmospheric medium. In the experiment, the effects of the 
ocean waves on vertical temperature profile and hence the refractive index profile is 
investigated. 

There is evidence from the experiment that the air turbulence due to the ocean waves 
resulted in deviations of the vertical temperature profile from the logarithmic profile thereby 
affecting the refraction of the optical signals. To analyze the refraction of the optical signal, 
only samples of those above the ocean waves, i.e., those not occluded by ocean waves, were 
used. For each time sample, the corresponding transmitter and receiver heights above mean 
sea levels were calculated from GPS surveying data and tide measurements provided by a 
National Oceanographic and Atmospheric Administration tide gauge. With the measured 
transmitter and receiver heights, and using the standard logarithmic temperature profile in the 
ray tracing program, it was observed that optical continuity between the transmitter and the 
receiver could not be established. An example of the ray tracing plot is shown in Figure 23. 
The mechanism in which the ocean waves affect the index of refraction needs to be examined. 
This research attempts to relate the effects of ocean waves on the temperature profile using 
the mixing length concept. 

The mixing length concept was introduced by Prandtl in 1925 to describe the 
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properties of turbulent transfer. In this model, it is assumed that there exists a length scale 
called the mixing length over which a turbulent eddy can travel and still maintain its identity. 
Details of the mixing length concept is contained in Appendix B. 

The effects of ocean surface waves on the gradients of the vertical temperature profile 
are discussed in the next section. 
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Figure 23. Optical Continuity cannot be established between Transmitter and Receiver. 



Source ht= 5.16 m 
Receiver ht = 4.69 m 
Range = 12.8 km 
Wave ht = 0.72 m 
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G. EFFECTS OF OCEAN SURFACE WAVES ON THE INDEX OF 
REFRACTION 

The refractive index at optical frequencies can be approximated by 



n - 1 =77.6(1 +7.52><10‘ 3 A" 2 )(P/J)10' 6 , 



(70) 



where P is the atmospheric pressure in millibars, T is the temperature in Kelvin and X is the 
wavelength of light in pm. In the atmospheric surface layer, pressure variation is small and 
can be neglected. Variation in the index of refraction is thus caused by the variation of 
temperature in the atmospheric surface layer. 

The vertical profiles of wind speed, u, temperature, T and humidity, q all have 
logarithmic forms £ =(u, T, q ), given by 





. \ Z K) 




(71) 



and the gradient profile given by 
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/ 



(72) 



th 

where is the flux scaling parameter, E, s is the value at the surface, is the roughness length, 
k is the von Karman constant and 4>^(z/Z)is the derivative version of the stability functions 
given in Appendix B. In Equation (72), the term kz is the turbulence mixing length. 

The above profiles are derived for turbulence above an immovable surface. A 
distinctive feature of the boundary layer adjacent to the ocean surface is the presence of a 
continuous random movement of the ocean surface. The presence of the random varying air- 
wave interface produces a constant variation in resistance of the ocean surface. The 
roughness elements that originate from the various components of the wave spectrum 
contribute in varying degree to the total wave dispersion. The air flow above the ocean waves 
induces a pressure field which transfers part of the energy of the air flow to the waves during 
wave development or conversely, the waves transfer energy to the air during swell 
attenuation. The surface waves, on the other hand, disturb the air flow above them and form 
air waves which interact with the gradient of the vertical profile. 

The wavy motion of the surface affects the profile deviation in two ways. Firstly, the 
short waves (high frequency component) contribute to the roughness of the water surface. 
The quantity z ^ specifies the scale of the turbulent eddies that are generated by the surface 
roughness elements. Experimental evidence [Ref 16] suggested that surface roughness should 
be proportional to wave steepness, HI A, where H is the wave height and A is the wavelength 
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of the dominant ocean wind waves, i.e., 






H K 
A g 



(73) 



As the ocean wave phase speed, C « [gA/(27t)] 1/2 , can be expressed as 



7 = a 

^0 u 



H_ 

2n 



\ 

K 

Cj 



(74) 



where a is the constant of proportionality. Equation (74) implies z 0u « H and z 0u « u t . A 
large u , implies a small T, as 



w'T' 






(75) 



and hence a smaller temperature profile gradient. In Equation (75), the overbar refers to 
averaging over a suitable time period. 

In this experiment, however, the wind speed was relatively low, with average wind 
speed of about 3.95 m/s. The dependence of roughness on wave steepness may be reduced. 
It may sufiice to use Chamock’s formula that has been interpolated to cover regime between 
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low wind and high wind speed [Ref. 16], 



z = 0.0185 — , 
e g 



(76) 



where g is the gravitational acceleration and u. is the friction velocity. The roughness length 
for a smooth surface is given by 




(77) 



where v is the dynamic viscosity of air. The total roughness length is given by 



Z 0u 



= z+ z. 



(78) 



The roughness lengths z 0T and z 0q are given by 



Z 0T 



0.4— 



u 

* 



(79) 
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and 



0 q 



v 

= 0.62— 



u 



( 80 ) 



Secondly, the longer surface waves do not contribute directly to the roughness, but 
they distort the air flow immediately above the ocean surface and are responsible for the 
vertical displacement of the frame of reference. Furthermore, the presence of the longer 
waves causes a redistribution of wave energy and the slope density of the short waves. This 
in turn influences the direction of the mean stress vector by creating an anisotropic roughness 
field at the surface. The long waves strained the shorter waves, causing them to turn in the 
direction of the longer waves. Donelan [Ref. 17] showed that the presence of long, gentle 
swell propagating in the wind direction has a significant effect on the wind sea spectrum. 
There is also evidence that counter swell, propagating against the wind led to higher drag 
coefficients than those experienced in pure wind waves [Ref. 14], Therefore the presence of 
swell will influence the turbulence of the air above it. By direct interaction with the air flow 
and due to the fact that swells travel faster than the wind, the swells will transfer momentum 
to the air flow above it and thereby enhance turbulent mixing. Since the random motion of the 
surface waves generate turbulence in the air above the surface waves and increase turbulent 
mixing, it is appropriate to employ a statistical mixing length scale that accounts for the 
increased mixing. The solution adopted here is to characterize the increased turbulence by an 
effective random mixing length scale that varies with the significant surface wave height, H s . 



70 



The significant wave height is one of the common parameter that oceanographers use as a 
measure for the average height of ocean surface waves. It is defined as the average of the 
highest one-third waves. The vertical profile of £ becomes 






In 



z+H. 



) 




( 81 ) 



and 



— = - 

dz k (z+H s ) \L 



( 82 ) 



Equation (82) indicates that the presence of ocean surface waves will reduce the profile 
gradient in a random manner that is determined by the surface wave height. An increase in 
wave height will lead to a reduction in the profile gradient. 

With the modified temperature profile, it was demonstrated by the ray tracing program 
that optical continuity between the transmitter and the receiver can be established. This is 
shown in Figure 24. This approach is based on applying the presumed scaling for the near 
surface condition. 
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Ray tracing using modified temperature profile, 30 May 96, 1225 hr 
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vn. CONCLUSIONS AND RECOMMENDATIONS 



This research showed that for optical signal samples that are not occluded by ocean 
surface waves, the signal fluctuations exhibit a log-normal distribution. This is characteristic 
of signals propagating through a random inhomogeneous medium. For signal samples that are 
occasionally obscured by ocean surface waves, the signal fluctuations deviate from a log- 
normal distribution. 

The effects of ocean surface waves on optical signal fluctuations are revealed in the 
temperature fluctuation spectral densities. The temperature spectral densities display 
prominent peaks at the corresponding peak frequencies of the ocean surface waves. These 
characteristic peaks are not present in the temperature spectral densities computed from 
meteorological measurements over land. 

Optical signal fluctuations increase as the signal ray path gets closer to the ocean 
surface waves demonstrating the significant effects of the ocean surface waves on signal 
fluctuations. Signal fluctuations were also shown to be dependent on wave age and the 
influence of ocean swells. Younger waves are rougher and are likely to cause a greater 
increase in signal fluctuations than more mature waves cause. When heights of signal samples 
are sufficiently low so that the ocean waves occasionally obstruct the ray paths, the ocean 
waves modulate the optical signals. Under moderate or high wind conditions, the ocean wind 
wave steepness should be incorporated into the roughness length scale, z Qu . 

On larger scale variations, the ocean surface waves modify the vertical temperature 
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profile and hence the refractive index profile. The standard logarithmic temperature profile 
derived for meteorological conditions over land is not suitable for the marine environment 
above the undulating ocean surface. Ray tracing using the standard logarithmic temperature 
profile indicates that optical continuity could not be established between the transmitter and 
the receiver. Unlike conditions over land, the ocean waves distort the air flow above and 
increase turbulent mixing. This increased mixing is proportional to the ocean wave height. 
The increased turbulent mixing can be accounted for by incorporating a statistical mixing 
length proportional to the ocean wave height. The proposed temperature profile assures 
optical continuity between the transmitter and the receiver. 

The research has demonstrated the significant effects of ocean surface waves on the 
propagation of low altitude optical signals above the ocean. The results of this study provide 
an insight into the effects of ocean surface waves on an optical signal skimming the ocean 
surface from an engineering view point. The study is not intended to disprove existing 
theories and hypotheses but rather to emphasize the relative importance of certain parameters 
that affect the propagation of optical signals in the vicinity of the ocean surface where the 
influence of ocean surface waves is significant. 

The existing model (e.g. ER. Tool) used by the Navy utilizes the standard logarithmic 
temperature profile to predict refraction of signals. The standard logarithmic profile is only 
valid for heights that are several significant wave heights above the ocean surface waves 
where the ocean wave effects are negligible. Within a few meters of the ocean surface waves, 
excessive ray bending occurs due to the logarithmic singularity of the temperature profile as 



74 



the ray approaches the ocean surface. The proposed statistical mixing length approach to 
modifying the temperature profile is perceived as an engineering fix to alleviate the large 
gradient in this critical layer above the ocean surface where sea-skimming threats are 
predominant. Further refinement in the modeling is recommended through extensive 
experimental measurements and data analysis. Modeling of signal propagation above an air- 
ocean interface spans an inter-disciplinary study of electronic engineering, oceanography and 
meteorology. The complexities of the air-ocean interface presents a unique challenge to 
scientists and engineers in a continuous search for more accurate modeling of signal 
propagation in this critical marine layer. 
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APPENDIX A. OPTICAL SIGNAL SAMPLES 
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Optical Signal Samples 





A 


B 


C 


D 


E 


F 


G 


H 


1 


1 


Date time 


Rx ht [m] 


Tx ht [m] 


height of 


rms 


% 


h A 2 


log var 


log var 


2 








ray path 


wave ht 








(wave occ) 


3 


5301225 


4.69 


5.16 


1.71 


0.36 


0.000106 


18.7128 


0.079 




4 


5301230 


4.56 


5.18 


1.64 


0.36 


0.000248 


16.5491 


0.082 




5 


5301233 


4.41 


5.19 


1.57 


0.36 


0.00062 


14.9584 


0.095 




6 


5301236 


4.27 


5.20 


1.50 


0.39 


0.005863 


14.5668 


0.067 




7 


5301838 


4.44 


4.94 


1.47 


0.7 


1.789256 


10.6063 




0.229 


8 


5301843 


4.26 


4.91 


1.36 


0.58 


0.956063 


7.9751 




0.198 


9 
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3.86 


4.82 


1.11 


0.45 


0.688336 


5.7584 




0.198 


10 
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3.77 


4.80 


1.05 


0.46 


1.123754 


5.1516 




0.145 


11 


5301909 


3.59 


4.78 


0.94 


0.63 


6.725038 


4.8771 




’ 0.089 


12 


5310943 


4.40 


4.95 


1.45 


0.2 


1.78E-11 


44.6558 


0.035 




13 


5310946 


4.24 


4.94 


1.37 


0.19 


3.23E-1 1 


39.5961 


0.065 




14 
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4.15 


4.93 


1.31 


0.22 


1 .22E-07 


36.3271 


0.123 




15 
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4.17 


4.87 


1.30 


0.25 


1.06E-05 


40.2262 


0.108 




16 
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4.10 


4.88 


1.27 


0.23 


1.88E-06 


37.5196 


0.143 




17 
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4.03 


4.88 


1.23 


0.2 


4.1E-08 


35.2933 
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18 


5311830 


3.86 


4.87 


1.13 


0.33 


0.030874 


15.0818 


0.154 




19 
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3.69 


4.84 


1.03 


0.32 


0.06726 


11.8792 


0.158 




20 
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3.53 


4.83 


0.93 


0.26 


0.016937 


9.2211 


0.157 




21 
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3.43 


4.82 


0.87 


0.34 


0.50895 


8.2193 


0.149 




22 
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2.94 


4.55 


0.48 


0.26 


3.26005 


3.2021 




0.087 


23 


5311937 


3.07 


4.53 


0.54 


0.29 


3.116286 


3.3635 




0.294 


24 


5311948 


3.32 


4.47 


0.66 


0.26 


0.586267 


4.4495 




0.244 


25 


6012208 


3.22 


4.05 


0.41 


0.26 


5.780919 


3.5517 




0.227 


26 


6012243 


3.21 


4.03 


0.39 


0.2 


2.511705 


3.3162 




0.252 


27 


6012308 


3.53 


4.05 


0.57 


0.19 


0.134146 


6.4756 


0.134 




28 


6021245 


4.39 


3.42 


0.57 


0.19 


0.010306 


12.4613 


0.257 




29 
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2.43 


3.53 


-0.23 
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85.2618 


3.0506 

7.2935 




0.086 


30 
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1.84 


3.69 


-0.52 


0.2 


99.52231 




0.103 


31 


6030930 


4.01 


4.01 


0.79 


0.22 


0.015119 


9.949 


0.107 




32 


6030937 


3.81 


3.96 


0.67 


0.22 


0.117065 


6.8136 


0.126 




33 


6031540 


4.26 


5.19 


1.49 


0.38 


0.004213 


24.3964 


0.103 




34 


6031551 


3.83 


5.23 


1.28 


0.3 


0.001032 


17.7704 


0.091 




35 


6031605 


3.34 


5.26 


1.02 


0.31 


0.052834 


11.5735 




0.132 


36 


6031611 


3.43 


5.28 


1.07 


0.34 


0.079321 


13.1283 




0.079 


37 


6041325 


3.58 


4.74 


0.92 


0.26 


0.020864 


11.8525 


0.157 




38 


6041417 


2.77 


4.68 


0.44 


0.34 


9.854486 


3.1197 




0.269 


39 
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3.38 


4.69 


0.79 


0.28 


0.249417 


7.3583 


0.145 




40 
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4.63 


3.38 


0.76 


0.31 


0.703241 


10.8039 


0.181 




41 
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3.42 


3.40 


0.20 


0.24 


20.57304 


1.5127 




0.196 


42 


6052000 


1.69 


3.42 


-0.72 


0.19 


99.9921 


8.4996 




0.123 


43 


6052041 


2.42 


1.69 


-1.17 


0.22 


99.99999 


19.4317 




0.077 
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Legend: 



Column A: 
Column B: 
Column C: 
Column D: 
Column E: 
Column F: 
Column G: 

Column H: 
Column I: 



Date and time of signal sample. 

Receiver height above mean sea level (in meters). 

Transmitter height above mean sea level (in meters). 

Lowest height of the ray path above mean sea level (in meters). 

Rms ocean wave height (in meters). 

Percentage of ocean surface waves protruding into optical signal ray path. 
Normalized mean square distance of optical signal ray path above ocean 
surface waves. 

Log-normal variance of signal samples not occluded by ocean surface waves. 
Log-normal variance of signal samples occluded by ocean surface waves. This 
is not strictly correct as the signal samples no longer exhibit log-normal 
distribution. 
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APPENDIX B. MIXING LENGTH HYPOTHESIS 



The mixing length concept was introduced by Prandtl in 1925 to describe the 
properties of turbulent transfer. In this model, it is assumed that there exists a length scale 
called the mixing length over which a turbulent eddy can travel and still maintain its identity. 

Suppose, under neutral atmospheric conditions, a parcel of air moves vertically 
through a small distance, /, from its original position, z+l, where the characteristic quantity, 
q(z), has an average value of q(z+l), then 

q(z) = q{z+t) . (83) 

Using Reynold’s decomposition, q(z) can be written as a sum of an average component and 
a fluctuating component, 

q(z) = q(z) + q'(z) . (84) 



Using mixing length hypothesis, 

?'<*) = «&♦/> - m - /I 2 . 

oz 



(85) 



Suppose — > 0, then as shown in Figure 25, / > 0 implies that the eddy comes from 

dz 

above and hence the vertical turbulent velocity, w' < 0 . By the same token, / < 0 implies 
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w' > 0. The above argument indicates that / and w' are negatively correlated, i.e., 



w'l < 0 . 



( 86 ) 



Assuming that mixing length hypothesis applies to the turbulent velocity component 
and the air temperature then 



T\z) 




(87) 



and 



// \ 1 du 

» 00 = 1 — , 

dz 



( 88 ) 



where T' is the turbulent air temperature and u' is the turbulent streamline wind velocity. 
Equation (88) implies 

— — —n du 

w u = w l — . (89) 

dz 

0 " 

In the atmospheric surface layer, — >0 and since w'l < 0, then 

dz 
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w'u' < 0 . 



( 90 ) 



Assuming the eddy velocity components, u' and w' have the same order of magnitude, then, 



/ / jdu 

w~ u ~ l — . 

dz 



Equations (89), (90) and (91) lead to 



/ 

w'u' = - / 2 

V 




2 



( 91 ) 



( 92 ) 



In the atmospheric surface layer, the mixing length varies with height above the earth’s 
surface. Near the earth’s surface, mixing cannot take place over large distances due to the 
solid boundary. Away from the surface, the turbulence is less inhibited and mixing is 
enhanced. In general, it is assumed that 

l = kz , ( 93 ) 



where k is the von Karman constant. 

In the atmospheric surface layer, the main task is to establish relationships between 
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vertical gradient of mean quantities of wind speed, temperature and specific humidity and 
surface fluxes. In particular, the surface fluxes of momentum, t, sensible heat, Q H , and latent 
heat, Q & are given by 



■ 2 

t = -p u'w' = pu, , 



(94) 



Qh = P C P T ' W ' = -P 



(95) 



and 



Qe = = ~P L e u *<l* ■ 



(96) 



where c p is the specific heat of air, and L e is the latent heat of vaporization of water. 
Combining Equations (92), (93) and (94), 



du = 
dz K z 



(97) 



where u is the mean wind velocity (the overbar has been dropped for convenience). Similar 
equations for T and q can be written: 
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8T = T\_ 
dz K z 



( 98 ) 



and 



dq = <L_ 
dz KZ 



(99) 



For non-neutral atmospheric conditions. Equations (97), (98) and (99) become 



du 

dz 



u. 



= — § u {z!L) , 

KZ 



( 100 ) 



dT 

dz 




( 101 ) 



and 
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( 102 ) 



I - - ■ 



KZ 



where (j>’s are the stability functions determined empirically, 



(1 -20 z/Ly m , 
1 +8 z/L 



z/L<0 
z/L > 0 ’ 



( 103 ) 



<i>r = ^ = 



| (1 - 16z/L)' 1/2 , 
[1+7 zIL, 



z/L< 0 
z/L> 0 



( 104 ) 



and z/L < 0, z/L >0, and z/L = 0 correspond to unstable, stable and neutral stratification 
respectively and L is the Obukhov length. 

Integrating equations (100), (101) and (102) yield profiles for mean u, T and q, 



2 / 



u = u s + — [\n (zlz 0u ) - iJj u (z/L)] , 

K 



( 105 ) 
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( 106 ) 



T=T s +^[ln(z/z 0T )-$ T (z/L)] 



and 



1C 



( 107 ) 



where u s , T s , q axe the values at the surface and i(r’s are the empirically derived stability 
correction functions given in Chapter V. 



87 







z + l 





z 


/ 
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q(z) 
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— Y 
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Figure 25. Mixing Length 
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m. OPTICAL SIGNAL SPECTRUM AND THE CORRESPONDING OCEAN 
SURFACE ELEVATION SPECTRUM 
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IV. RAY TRACING USING STANDARD LOGARITHMIC TEMPERATURE 
PROFILE 
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V. RAY TRACING USING MODIFIED TEMPERATURE PROFILE 
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